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Abstract. 

We present a novel laser-based ion acceleration scheme, where a petawatt 
circularly polarized laser pulse is shot on an ultra-thin (nano-scale) double-layer 
target. Our scheme allows the production of high-quality light ion beams with 
both energy and angular dispersion controllable by the target properties. We show 
that extraction of all electrons from the target by radiation pressure can lead to a 
very effective two-step acceleration process for light ions if the target is correctly 
designed. Relativistic protons are predicted with pulse powers of a few petawatt. 
Careful analytical modeling yields estimates for characteristic beam parameters 
and requirements on the laser pulse quality, in excellent agreement with one and 
two-dimensional Particlc-in-Cell simulations. 
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1. Introduction 

Generation of high-energy ion beams by interaction of an ultra intense laser pulse with 
a solid target is one of today's hot topics in laser-plasma interaction. Such ion beams 
have properties making them very interesting for a wide range of applications, such 
as proton radiography [l] , fast ignition in the context of inertial fusion [2 [3l |4] , or 
hadron-therapy [5] . While generation of ions with energies up to several tens of MeV 
has already been demonstrated [5], controlling their energy distribution remains a 
crucial issue for most applications. 

Different mechanisms of ion acceleration have been proposed depending on 
whether ions originate from the front-side (irradiated by the laser) or the rear-side of 
the target. Ion acceleration at the target front-side occurs mainly in the electric field 
resulting from electron sweeping at the front of the laser pulse [7] , leading to creation of 
collisionless electrostatic shocks [8l[9] or solitary waves [10]. At the target rear-side, ion 
acceleration occurs in the strong electrostatic field resulting from charge separation 
due to hot electrons escaping into vacuum. This mechanism, referred to as target 
normal sheath acceleration (TNSA) [11], is the dominant process of ion acceleration 
for currently available (moderately relativistic) laser intensities [12]. TNSA provides 
ion beams with interesting properties, such as good laminarity, small aperture angle (a 
few degree) and relatively large efficiency of energy conversion (a few percents) from 
the laser pulse to the ions. However, the resulting ion beams have a characteristic 
broad (quasi-thermal) spectrum with a sharp cut-off at maximal energy. 

Several proposals have been made how to control the energy distribution of 
laser-created ion beams. Recently, efficient acceleration with relatively small energy 
dispersion has been observed in numerical simulations where a circularly polarized 
(CP) laser pulse was focused on a thin target [131 El HSl HSl HI]- Using CP light 
indeed allows to strongly reduce electron heating [5] and therefore prevents TNSA 
and the associated broadening of the ion spectrum. Ion acceleration then follows from 
a front-side mechanism referred to as laser-piston or light-sail acceleration. The whole 
target is accelerated as a neutral bunch, resulting in a quasi-monochromatic ion energy 
distribution. However, due to the non-homogeneous field in the piston, the presence 
of low-energy ions in the beam cannot be avoided. 

Moreover, alternative methods based on multi-species (homogeneous or 
multilayered) targets have been proposed to control the ion spectrum [18l [191 (iQl 
[211 [22l [23] . In particular, double layer targets have attracted a lot of interest as they 
allow to increase ion beams monochromaticity |191 120j. In the proposed schemes, a 
laser pulse is focused on the first target layer that consists of heavy, highly charged, 
ions. Electrons gain energy in the laser field and are either heated or completely 
extracted from the laser focal spot. In the first case, ions of the second layer are 
accelerated in the ambipolar electrostatic field created at the rear-side of the first 
ion layer by hot electrons [21]. The possibility to generate 1.3 MeV proton beams 
with energy dispersion ~ 25 % and 3 MeV carbon beams with energy dispersion 
^17% using double layer target has already been demonstrated experimentally by 
Schwoerer et al. |19j and Hegclich et al. [20], respectively. The mechanisms behind 
these observations are similar to TNSA, and generation of hot electrons is a dominant 
process under present experimental conditions. This rather complex mechanism makes 
the control of the ion beam properties non trivial. In the second case, where electrons 
are removed from the target, light ions are accelerated in the strong electrostatic field 
of the expanding first ion layer. This is the so-called regime of directed Coulomb 
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explosion (DCE) [T51 [21]. While ion acceleration in the ambipolar field is highly 
sensitive to the hot electron temperature and thus to the laser parameters, ion 
acceleration by DCE depends only on the target properties, thus allowing a better 
control of the ion source. Nevertheless, no efficient mechanism has yet been proposed 
to expel all electrons from the target and achieve optimal DCE. 

In this paper, we present an elaborated laser-based ion acceleration scheme which 
allows to control energy dispersion, potentially below 10 %. This scheme relies on the 
complete, laser-induced, removal of electrons from an ultra-thin double layer target. 
To make electron sweeping by the laser radiation pressure more effective, we make 
use of a CP laser beam. One-dimensional (ID) and two-dimensional (2D) numerical 
simulations are performed with the particle-in-cell (PIC) code CALDER [25]. They 
show that complete extraction of electrons from the target is possible above a 
threshold intensity that depends mainly on the target areal charge. The resulting 
ion acceleration from the double layer target is then discussed using both numerical 
simulations and analytical modeling. First, light ions making the second target layer 
are accelerated in the quasi-homogeneous electrostatic field created between the first 
ion layer (that consists of heavy, highly charged, ions) and the forward-going electron 
cloud. This first stage, referred to as linear plasma acceleration (linPA), ends when 
electrons are pushed far enough from the heavy ion layer. The subsequent acceleration 
phase of light ions depends on whether they have acquired or not relativistic velocities 
in the linPA stage. Especially, it is shown that low-energy ions can gain further energy 
in a stage similar to DCE, while high-energy, relativistic, ions gain most of their energy 
in the linPA stage of acceleration. The final energy and energy dispersion are shown 
to depend mainly on the target composition and geometry. 

Because this novel mechanism of light ion acceleration depends mainly on 
the target properties, it is thus much less sensitive to the inevitable shot-to-shot 
fluctuations of the laser parameters than TNSA-based schemes. The ability to control 
the ion beam properties by a careful design of the target makes this acceleration 
process particularly interesting for applications requiring high-quality ion beams, such 
as hadron-therapy. 



2. The linear plasma accelerator (linPA): Theory 

Our characteristic target (see Fig. [T]) consists of a first layer of heavy, highly charged 
ions with mass mi 3> mp (rrip ~ 1836 is the proton mass) and charge Z/j ^ 1, atomic 
density and thickness dh- In this paper, mass, densities, charges and distances 
are normalized to the electron mass rrie, critical density ric = eoWeW^/e^ (eg is the 
permittivity of vacuum, e is the electron charge and cjl is the laser frequency), electron 
charge e and inverse laser wave number k^^ = c/lol, respectively. A second layer, with 
density ni and thickness di, contains the light ions whose acceleration is considered. 
As will be discussed later in the paper, these ions must have a charge-over-mass ratio 
larger than the one of the so-called heavy ions. This is usually the case for species 
with mi <C m/j, first because of the neutron contribution to the mass of the nucleus, 
and second because the heavy ions may not be fully ionized. This target is irradiated 
by a CP laser beam at relativistic intensity. In what follows, the laser field amplitud^ 
ql ^ El/ Ec > 1, as well as all other electric fields, are given in units of the Compton 

X Normalization arc chosen so that = 1.38 X 10^* W/cm^ Mm^, with the laser intensity 

and the laser wavelength. 
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Figure 1. Characteristic double-layer target used in this paper and spatial 
distribution of the electrostatic field along the x-direction (green curve). 
Quantities ui^i dh and uij, d; denote the transverse size and thickness of the 
heavy and light ion layers, respectively. The electron bunch (on the right) has 
the transverse size wj_ ~ min{ui(,, ui^}, where wi^ is the width of the laser focal 
spot. The laser pulse propagates from left to right. 
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2.1. Electron response to laser radiation pressure: creation of the linPA 

In a first stage, electrons are quickly pushed forward into the target by radiation 
pressure 11^ = (1 + i? — T) a\/2, where 11^ is given in units of rig rUe c^, and R and 
T are the laser reflectivity and transmittance. As electrons penetrate deeper into the 
target, an electrostatic field is built up at the target front side. If the target is thick 
enough, a quasi-equilibrium between electrostatic and radiation pressures can set in. 
This process has been observed in different studies related to hole-boring [SI [3] or 
hght-sail acceleration [13l [HI [IB Ull [H [22| , using both linearly polarized (LP) or CP 
hght. 

In this paper, we consider the case where radiation pressure is too high to be 
balanced by the maximum electrostatic pressure a^otl'^ (where Otot =■ a/i + o-i is the 
maximum electrostatic field due to complete electron-ion separation; = Zh dh 
and ai = Zi ni di are the electrostatic fields due to the bare ions of the first and second 
layers, respectively). When the laser field amplitude ol becomes larger than aatot, 
with a of the order of unitjH, electrons are completely extracted from the targelljl. 
To our knowledge, this behavior was reported for the first time in Ref. [31] where 
the possibility to generate dense electron bunches with ultra-intense laser pulses with 
sharp rising edge is discussed. More recently, Rcfs. [161 [2] have reported similar 
electron behavior when investigating the transition from RPA to Coulomb explosion 
of a monolayer target. It is interesting to note that, for a given target electron density, 

§ The parameter a on laser intensity accounts for relativistic corrections, spatial and temporal profile 
of the laser pulse, etc. 

II Let us note that, considering a plasma with ion density <C m/Z^ and a laser field amplitude <C m 
(where m and Z are the ion mass and charge, respectively), the ion layer is transparent to the laser 
beam. 
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complete extraction of all electrons indeed requires the target thickness to be smaller 
than the so-called optimal thickness for RPA. 

For laser field amplitude 3> 1, electrons become relativistic almost 

instantaneously and the separation time scale is tc ^ dh + di (in units of i^^^)- This 
time scale is much shorter than the pulse duration and all characteristic times of ion 
acceleration. 



2.2. Ion acceleration in the linPA 

Right after the electrons are separated from the ions, the target exhibits a capacitor- 
like structure, with a uniform electrostatic field of amplitude at^t built up between the 
electron and ion layers. In this Section, we restrict our study to a ID problem, thus 
the electrostatic field amplitude does not depend on the distance between the ion and 
electron layers. Limiting multi-dimensional effect are discussed later in Scc|4l 

The electric field seen by the heavy ions of the first layer is not homogeneous. 
Assuming a fiat-top ion density profile, it increases linearly from to so that the 
velocity of a given ion depends on its initial location, thus yielding a large energy 
spread. If the charge-over-mass ratio of the light ions is larger than the one of heavy 
ions, there are no intersections between the trajectories of heavy and light ions, and 
the electrostatic field in the second ion layer varies linearly from up to at^t. With 
O'l ^ o,h-: this field is quasi-homogeneous and small energy dispersion for light ions is 
expected. 

In ID geometry and assuming complete electron expulsion, light ions see a 
constant, quasi-homogeneous, accelerating field. Relativistic equations of motion for 
the slowest and fastest light ions, experiencing the accelerating fields ah and atot 
respectively, can thus be solved exactly. Extracting the ion mean energy £ and energy 
dispersion AS at time t after creation of the linPA is straightforward: 



£ ^ mi 



vi + ^W-i , (1) 

(2) 



^l + {l + ai/anft^/tl-^l + t^/tl 



where energies are normalized to the electron rest energy me and times to ijjj} ■ 

The characteristic time tr = mi/{Ziah) denotes the time required for light ions 
to gain relativistic energies. In the limit t <^ tr, light ions have non- relativistic 
velocities and the mean energy and relative energy dispersion are £^'^^ ~ m;<^/(2i^) 
and (Af )'^^^ ~ 2 ai/au, respectively. On the contrary, if the linPA can be maintained 
over times t ^ tr, light ions obtain ultra- relativistic velocities. Their mean energy then 
evolves as ~ mi t/tr, while their relative energy dispersion is (A^^/f)'"'') ~ ai/oh. 
In both limits, the relative energy dispersion remains small if the electrostatic field ai 
due to light ions is small compared to the accelerating field at (due to heavy ions), 
and it can be controlled by adjusting the target properties. The acceleration scheme 
then suffers limitations similar to any linear accelerator: controlling the relative energy 
dispersion limits the areal density of accelerated ions ctx = ni di (in units of nc/k^). 
Nevertheless, it is shown in what follows that satisfactory values of a± can be obtained 
using laser field amplitude ~ 100. 

It is interesting to note that the characteristic time tr as well depends only on 
the target parameters. In particular, for a given charge-over- mass ratio Zi/mi, tr 
depends only on the amplitude of the accelerating field ah. Obviously, the larger ah, 
the shorter tr and the higher the light ion energy at a given time t. Moreover, we 
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want to emphasize here that the main hmitation to the ion energy follows from multi- 
dimensional effects. While these effects are discussed in more detail in Sec. 21 it is 
worth pointing out at this stage that they set in after electrons are pushed on a distance 
larger than either the laser focal spot diameter or the transverse width Wh of the 
first target layer. Assuming relativistic electrons, this defines a time tunPA S "^i- 
(where w± = uii-a{wL,Wh}) during which light ions can be efficiently accelerated in 
the electrostatic field ah ■ The condition for obtaining relativistic ions on a time tunPA 
defines a minimum value for the accelerating field a^f^'' ~ mi/ {Ziw±_). Considering 
that the condition for charge separation (cll > cecLhl'^) defines a threshold for the laser 
field, and noting that = a| w'^ defines the effective laser power on the target (in 
units oi Ic /k\), one can define a characteristic laser power P^j//^ = mf /Zf required 
to obtain relativistic light ions in the linPA. This characteristic power depends only 
on the light ions charge-over-mass ratio and it is typically of a few PW for protons. 

3. The linPA: One-dimensional numerical simulations 

3.1. Numerical results vs. theoretical predictions 

One-dimcnsional simulations using the PIC code CALDER have been performed to 
explore this new regime of ion acceleration and to test theoretical predictions. In these 
calculations, a CP laser beam with field amplitude = fOO {II ^ 1.4 x 10^^ W/cm^ 
at a wavelength Ai, = 1 /im) is focused at normal incidence on an ultra-thin (few 
tens of nm) double layer target. In a first attempt, a flat-top temporal laser intensity 
profile is chosen in order to simplify comparison to analytical estimates. 

The first target layer is made of carbon (Z^ = 6, to/,, = 12 to^) with atomic density 
nil = 58 6.4 X 10^^ cm~^) and thickness ranging from = 0.04 up to 0.12 (6.4 to 
19.1 nm). The second layer contains only hydrogen (Zi = 1, mi = mp) with density 
m = 5.8 (~ 6.4 X 10^^ cm-3). Its thickness, di = 0.04 - 0.7 (6.4 to 111 nm), is adjusted 
to control the energy dispersion of the resulting ion beam as well as its areal density 
a± = nidi. The foil has step- like density profile in each layer. We assume complete 
ionization at the beginning of the simulation and the initial electron temperature is 
1 keV. 

In the following simulations, the numerical domain is 7 Xl long with mesh size 
dx = 2.5 X 10-3 ^ A^'^, where A^^ = \/Te/{Zh 71 /i) is the normalized Dcbye length in 
the first target layer and Tg ~ 1/511 is the normalized initial electron temperature. 
The simulation duration is 5 laser periods and the time step dt dx/2. The laser 
propagates from the left to the right and reaches the target ^ Itl after the beginning 
of the simulation (tl is the optical cycle). This time is referred to as the time 
zero. Entrant (absorbing) and absorbing (absorbing/reinjecting) boundary conditions 
are used for the electromagnetic field (particles) at the left and right edges of the 
simulation box, respectively. 

Figure [5] shows a snapshot of the ion and electron densities and the electrostatic 
field four laser periods after the beginning of the interaction, in the case where 
di = 0.235 37.4 nm) and for two different thicknesses of the first layer c?/, = 0.08 
(~12.7 nm) and dh = 0.12 ('^19.1 nm), see Figs. [2^1 and b, respectively. We can 
observe that, in the case of the thinnest target (Fig. [2^), electrons are pushed as a 
compact bunch. As confirmed in the phase-space in Fig. [2j;, all electrons are pushed 
forward at relativistic velocities. For the thickest target, however (Figs, [^b andl^Ji), 
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Figure 2. a) Snapshot at 4 = 4t£ after the beginning of the interaction of the ion 
(blue: carbon, red: hydrogen) and electron (black) density and electrostatic field 
(green). The laser propagates from left to right. The target thickness is dh = 0.08 
(~12.7 nm). b) Idem but with the target thickness = 0.12 (~19.1 nm). 
c) Electron momentum in the direction of the laser propagation versus position 
for parameters of panel a. d) Idem but for parameters of panel b. 



electrons are not pushed as a compact bunch anymore. Figure [2]i shows that, in 
this case, the electron dynamics is rather complex. Because the laser field does not 
strongly exceed the electrostatic field due to electron-ion separation, some electrons 
are pushed forward by the laser pulse while others are accelerated backward in the 
strong electrostatic field. A similar behavior has already been reported in the laser- 
piston regime by Naumova et al. [4]. As demonstrated by these authors, the backward 
accelerated electrons may interact with the incident laser pulse, experience important 
radiative friction |32] and thereafter be slowed down. However, this effect is not taken 
into account in our simulations. Later in this Section (see Sec. 13. 2p . more details are 
provided concerning the electron behavior and their effects on ion acceleration. 

Figures[2^ and[2]3 also show that carbon ions are accelerated in an inhomogencous 
(linearly varying) electrostatic field, and hence get smeared out along the laser 
propagation direction. On the contrary, protons see a quasi-homogeneous field, whose 
amplitude, ~ 28 (- 90 TV/m, see Fig. Hi) and - 42 (- 130 TV/m, see Fig.^j), is 
in excellent agreement with the analytical prediction atot o-h ~ dh- Moreover, 

because the charge-over- mass ratio is larger for protons than carbons, the proton layer 
is quickly separated from the carbon one. 

Figure [3] summarizes the proton beam properties obtained in numerical 
simulations and compares them to theoretical predictions. Temporal evolution of 
proton energy is shown in Fig. [3^ for different thicknesses of the first layer for 
parameters previously detailed. After a few laser cycles only, proton energies up 
to 150 McV are obtained, in excellent agreement with analytical predictions from 
Eq. ([1]). Complementary simulations are also presented considering a carbon foil with 
thickness dh = 0.27 43 nm) and solid density rih — 91 (10^^ cm~'^) irradiated by 
a CP laser with field amplitude a^, = 400 (orange curves). Under such conditions. 
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Figure 3. a) Temporal evolution of the ion energy for different first layer 
thicknesses, df^ = 0.04 (~6.4 nm) (green, triangles), df^ = 0.08 (~12.7 nm) (black, 
circles) and df^ = 0.12 (~19.1 nm) (blue, square) and for parameters allowing 
relativistic ion generation (orange, star). The right vertical axis provides the 
energy corresponding to this last simulation. Dashed curves account for analytical 
estimates from Eq. JTJ. b) Relative energy dispersion A£/£ in % versus the 
second layer thickness di for target parameter corresponding to panel a. c) Energy 
spectrum of protons in the case where df^ = 0.08 (~18.8 nm) and for three different 
second layer thicknesses: d; = 0.05 (~6.2 nm) (solid line), d; = 0.25 (~30 nm) 
(dashed line) and di = 0.5 (~60 nm) (dotted line), d) Dependence of the areal 
density of the proton beam on the second layer thickness d; . 



protons undergo a strong acceleration and can gain relativistic energies in only a few 
optical cycles. Here, protons with energy larger than 1 GcV arc obtained, as expected 
from analytical modeling. This confirms that the linPA allows to generate high energy 
ion beams with high power laser on only a few laser cycles, that is without requiring 
exceedingly large laser energies. 

As previously underlined, the energy dispersion of the ion beam can be tuned by 
controlling the ratio ai/ah, i.e., the target properties. Figure [3)d presents the energy 
dispersion obtained in numerical simulations for different ratios ai/ah- As expected, 
the proton beam monochromaticity can be considerably improved by decreasing the 
thickness of the second layer, whereas the ion energy, which depends mainly on the 
first layer properties, is not modified. Energy dispersions of the order of 2 % are 
obtained, which is very attractive for possible medical applications. Moreover, we 
point out that estimates from Eqs. ([T]) and ^ compare well with numerical results. 
Actually, they may even overestimate them in those cases where some electrons partly 
neutralize the proton bunch, reducing the field inhomogeneity (Fig. Os). Last but not 
least, we want to emphasize that, in contrast to acceleration by radiation pressure (see 
e.g. Refs. [m [m Hi]), and as can be observed in Fig. [3h, all protons are contained 
within the monochromatic peak. 
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The drawback of the low energy dispersion permitted by this method is the 
low areal number of accelerated protons. However, as shown in Fig. for typical 
parameters of this study, a± = nidi up to a few 10^ particles//im^ can be reached, 
while keeping energy dispersion to a few percent only. 

To demonstrate acceleration of light ions other than protons, ID simulations 
have been performed to investigate the generation of quasi-monochromatic carbon ion 
beams. Let us consider a Al = 0.52 /im-laser pulse with field amplitude aL = 300 
{II 3.3 X 10^2 W/cm^) focused onto a 5 nm-thick {dh = 0.06) gold target at solid 
density rih = 14.8 (~ 5.9 x lO^^cm""^). Because ionization is not accounted for in 
our code, we discuss two sets of simulations where gold is either completely ionized 
{Zfi = 79 leading to ah ~ 70.2) or electrons of the three inner shells remain bound 
{Zh = 51 leading to ah = 45.3). A fully ionized, thin, carbon layer with density n/ = 25 
(~ lO^'^cm"'^) is placed at the rear side of the gold layer. Its thickness is adjusted so 
that relative energy dispersion of the carbon ion beam remains ^ 10%: di ~ 0.045 
(~ 4 nm) for the Au™+ layer, and di = 0.030 (~ 2.5 nm) for the Au^^+ layer. As 
a result, the generation of carbon beams with relative energy dispersion of the order 
of 15%, and energy 190 MeV/nucleons (using Au''^"'") and 80 MeV/nucleons (using 
Au^^+), is observed, in good agreement with theoretical predictions from the linPA 
model. This confirms that the proposed mechanism of acceleration can be applied to 
ions heavier than protons. However, it requires that the ions of the first target layer 
have a smaller charge-over-mass ratio than the ions whose acceleration is considered. 
Here, the gold ion charge-over-mass ratio is ~ 0.40 and ~ 0.26 for Au''^+ and Au^^+, 
respectively, and 0.5 for the fully ionized carbons. 

3.2. Influence of the laser temporal profile and polarization 

So far, simulations have been performed using an instantaneous ramp-up of the laser 
intensity. However, electron response to the laser pulse may be strongly dependent 
on its temporal profile. We have therefore tested the effect of a finite rise-time of the 
laser intensity on the electron dynamics and the subsequent light ion acceleration. For 
this purpose, a target with thicknesses dh = 0.08 (12.7 nm) and di ~ 0.16 (25.4 nm) 
irradiated by varying rise-time laser pulses is considered. Figure 3] shows the temporal 
evolution of electron density along the laser propagation axis x. For times i < 0, the 
laser pulse has not yet reached the target, which undergoes expansion due to non-zero 
initial temperature only. For times t > 0, electrons are strongly accelerated forward 
by the laser pulse. Once they leave the target, the electrostatic field atot is built 
between the electron cloud and ion layers. For sufficiently long times, this field may 
be strong enough to accelerate backward part or most of the electrons, as previously 
observed in Fig.[2]3. Considering a laser field amplitude = 100 and a flat-top laser 
intensity profile (Fig. |4^), a small fraction of the electrons returns toward the ion 
layers 5tl after the beginning of interaction. This time is reduced when the linear 
ramp-time is increased. For a 2tl ramp-time and laser amplitude ol = 100 (Fig. 0)3), 
a considerable fraction (> 40 %) of the electrons is accelerated backward less than 4 tl 
after the beginning of interaction. Increasing the ramp-time up to 4 tl even makes it 
impossible to extract all electrons from the ion layers, thus preventing efficient light 
ion acceleration in the linPA. 

Therefore, our proposed mechanism for ion acceleration is very sensitive to the 
laser pulse profile. It is nonetheless possible to mitigate these effects by increasing 
the peak laser intensity so as to steepen the intensity temporal profile. Figures |4ji 
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Figure 4. Contour plot of the electron density along the laser propagation axis 
X and time. The CP laser pulse has field amplitude = 100 (a,b,c), = 200 
(d,e,f) and = 400 (g,h,i). Without ramp time {a,d,g), with ramp time 2tl 
(b,e,h) and with ramp time 4t]^ (c,f,i). 



and [3^ indeed show that eleetrons ean be pushed forward as a eompressed buneh 
during the whole simulations by increasing the laser field amplitude up to = 200 
and aL = 400, respectively. More precisely, for a linear ramp-time 2tl, only a few 
electrons are accelerated backward when ul = 200 (Fig. while all electrons are 
pushed forward for ol = 400 (Fig. [4ji). However, an efficient piston- like acceleration 
remains difficult to achieve for a 4 ramp-time (Fig. |4f ) since it requires a laser field 
amplitude aL — 400 to prevent electrons from being accelerated backward (Fig. |4l). 

The influence of the laser polarization has also been investigated. The temporal 
evolution of the electron density along the laser propagation axis, x, when irradiated 
by a LP laser pulse is presented in Fig. [5] for field amplitudes ul = 100 and ul = 400, 
and for different ramp-times. When irradiated by LP light, electrons are strongly 
heated [9l[T4l[27j. To remove all electrons from the target, radiation pressure has to 
overcome not only electrostatic pressure but also thermal pressure due to extremely 
hot electrons. Much higher laser intensities are therefore required, which explains why 
the linPA regime was not relevant under conditions of Ref. [18]. At = 100 (Figs. [5^), 
some electrons are accelerated backward after only ~ 2.5 tl, as compared to ~ St/, 
when considering a CP laser pulse with similar parameters (Fig. [4^). Enhancing the 
linear ramp-time (Figs.[5}3 and[5j:) further hinder the electrons' removal by enhancing 
thermal pressure. This leads us to the conclusion that, for parameters characteristic 
of this study, and for a given laser intensity, using CP light remains more efficient to 
extract electrons from the target. 

As can be observed in Figs. |4] and EJ electrons display a very complex dynamical 
behavior. In particular, the early stage of the laser pulse interaction with the target 
is characterized by a strong acceleration of electrons, which gain relativistic energies 
in less than one optical cycle. Under such conditions, the Doppler-shift is modified 
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Figure 5. Contour plot of the electron density distribution along the laser 
propagation axis x and time. The LP laser pulse has field amplitude = 100 
(a,b,c) and = 400 (d,e,f). Without ramp time {a,d), with ramp time 2 ti^ (b,e) 
and with ramp time At]^ (c,f). 

and therefore the radiation pressure. Furthermore, due to the large field amplitude 
(ol > 100 typically) and the small target thickness, the electron cloud is not opaque 
to the laser. Relativistically induced transparency occurs |28| . and the laser pulse 
penetrates through the whole electron cloud, even though this cloud expands on 
several wavelengths. Obviously, this also modifies the radiation pressure. While these 
phenomena are accounted for in our PIC simulations, radiation friction is not taken 
into account. It may also strongly influence the electron dynamics and in turn ion 
acceleration, as previously observed in Ref . [D . Modeling the electron dynamics under 
such conditions thus appears to be very challenging and we leave this complex issue 
for future investigations. 

To conclude with the electrons' behavior, it must be stressed that both ID and 
2D simulations overestimate the electrostatic field between the electron and ion layers, 
thus causing the electrons to remain bound to the ion layers. Conversely, in 3D 
geometry, electrons are able to separate definitely from the ion layers when their 
kinetic energy Uk = "fe ~ ^ becomes larger then their potential energy in the ions field. 
Up = Q/{47Td), where Q ^ ahw\ is the total charge of the electron cloud in units of 
e nc/k\ and d is the distance between the electron cloud and the heavy ion layer |29j . 

From previous simulations, one can expect electrons to reach kinetic energies 
10 < C/fe < 100 ('^ 5 — 50 MeV) after only a few optical cycles in the laser field (e.g. 
see Fig. [2}:). Considering that electrons can be efficiently pushed over distances of a 
few wavelengths on such time scales, generation of electron bunches with the charge 
a few nC is expected under current conditions of irradiation and target thickness. 

4. Multi-dimensional effects on the linPA and directed Coulomb 
explosion (DCE) of the target 

Previous modeling and numerical simulations rely on an idealized ID picture of the 
linPA. In this Section, 2D simulations are discussed that show how multi-dimensional 
effects arise when electrons are pushed far enough from the ion layers. Their impacts 
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on ion acceleration are discussed. Especially, we show that linPA is only the first stage 
of the acceleration process. A second acceleration stage occurs when the distance 
between electrons and ions becomes similar to the transverse size of the system. 

Two-dimensional numerical simulations 

Two-dimensional simulations have been performed with CALDER to investigate multi- 
dimensional effects on the proposed acceleration scheme. In these simulations, the 
CP laser pulse has a flat-top temporal profile with duration Tp = 10 (~ 33 fs), 
and the maximum field amplitude is = 100 (~ 1.4 x lO^^W/cm^). Along 
the transverse j/-direction. the laser intensity follows a 6th-ordcr super-Gaussian 
distribution. Two value of the laser focal spot full-width at half-maximum have been 
considered: wl = 10 and wl — 20 Al 10 — 20 /im, respectively). This laser pulse 
is focused on an ultra-thin double layer target following the design considerations of 
the previous sections. The target's first layer is made of carbon with atomic density 
n/i = 58 6.4 X 10^^ cm^^), thickness dh ~ 0.08 (~ 12.7 nm) and transverse width 
Wh = 200 30 /im). An hydrogen dot with thickness di = 0.16 25.4 nm) and 
transverse width wi = A^ (^^ 1 /im) is placed on the rear-side of the first layer. Two 
different atomic densities, ni = 6 and n; = 12 0.6 — 1.2 x 10^^ cm^"^, respectively), 
have been considered. 

The simulations results are given in Figs. [S] and \7\ Figure [S] shows the spatial 
distributions of electrons, carbon ions and protons at three different times with the 
parameters ni ~ 6 and wl = 10 A^. In Fig. [6^, the characteristic features of the 
linPA arc clearly apparent. Four laser periods after the beginning of the interaction, 
an electron bunch with the transverse width w± ^ 12 A^, of the order of the laser 
focal spot, is separated from the ion layers. It is pushed forward over a distance of the 
order of 4 Al (consistent with the previous consideration that electrons leave the target 
with the light velocity). A strong electrostatic field, at ~ 28 90TV/m), is built 
up between the electron and carbon layers. This is confirmed in Fig. [7^, where the 
longitudinal component of the electrostatic field at the position of the proton bunch is 
shown as a function of time. The early stage of proton acceleration is characterized by a 
constant accelerating field, and its duration, tunPA, can be extracted from Fig. [7^. For 
a laser pulse transverse width — 10 Al, one has tunPA ~ 7 A^, while tunPA ^ 9.5 \l 
for Wl — 20 Al. At this point, it is interesting to note that tunPA is not twice as long 
for WL = 20 Al as for wl = 10 Al. This is because in the case where wl = 20 Al, 
some electrons come back toward the ion layers and thus reduce the accelerating field 
prematurely. For short times, t < tunPA, protons are consequently accelerated in 
the capacitor-like electrostatic field characteristic of the linPA process. Figure [TJd 
shows the temporal evolution of the proton energy for different target parameters. As 
predicted from the linPA model, the proton energy during this stage depends neither 
on the second layer's density ni, nor on the first layer's width. In this early stage, an 
excellent agreement with theoretical predictions is obtained. 

In addition. Figs. ^ and [TJi show the temporal evolution of the proton bunch 
energy and angular dispersions. Especially, Fig. [Tj; shows that, while the proton 
energy dispersion is not strongly modified by increasing the laser pulse transverse 
width, it can be tuned by changing the density of the hydrogen dot. This confirms 
that the energy dispersion follows from electrostatic repulsion inside the proton bunch 
itself. Energy dispersions of the order of 7 % and 14 % are obtained depending on the 
density ni, in very good agreement with predictions from the linPA model. Similar 
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Figure 6. Snapshots of the electron distribution (gray scale), carbon contour 
plot (blue) and proton contour plot (red) for different times: a) t = 4r^ after the 
beginning of the interaction; b) t = 14 tl after the beginning of the interaction 
and c) t = 30 after the beginning of the interaction. 



observations follow from Fig.[7ti concerning the angular dispersion of the proton bunch. 
In the linPA stage of proton acceleration, the angular aperture of the proton beam is 
therefore mainly governed by electrostatic self-repulsion. For parameters of this study, 
angular aperture of a few degree only are obtained. 

Let us now investigate proton acceleration on times larger than tunPA- As can 
be observed in Figs.[6lD and[6t, electrons are here pushed over a distance of the order 
of, or larger than w± . Multi-dimensional effects thus set in, reducing the accelerating 
field (Fig. [7^) and making proton acceleration less efficient (Fig.[7)D). While the proton 
energy in the linPA stage of ion acceleration did not depend on the transverse width wl 
of the laser pulse, increasing allows to enhance the final proton energy by delaying 
multi-dimensional effects. Furthermore, in this second stage of proton acceleration, 
one can observe an enhancement of both the energy and angular dispersion of the 
proton bunch (Figs.[71: and[71i). While energy dispersion is still mainly determined by 
the charge in the second layer (Fig. [71:), one can observe that the proton beam angular 
aperture depends on both the carbon layer's width and the density of the hydrogen dot 
(Fig. [31). This prompts us to suggest that, in this phase of the acceleration process, 
both the transverse inhomogeneity of the accelerating field and Coulomb self-repulsion 
are responsible for the proton bunch angular aperture. At the end of the simulation, 
energy dispersion of the order of 10 % and angular aperture < 3° are obtained. 

These simulations demonstrate that even when multi-dimensional effects are 
accounted for, the linPA mechanism is effective for accelerating light ion beams, as the 
first, transient, stage in a two-step acceleration process. In 3D geometry, the duration 
of the linPA stage might be further reduced thus limiting the final proton energy. On 
the other hand, one expects 3D effects to mitigate the energy and angular dispersions. 

4-2. Estimates for three-dimensional effects 

To understand the behavior of light ions over times larger than tunPA , it is necessary 
to compute the distance traveled during the linPA stage. The position of the slowest 
and fastest light ions at t = tunPA can be easily obtained by integrating the relativistic 
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Figure 7. Temporal evolution of: a) the accelerating electrostatic field along the 
x-direction at the position of the proton bunch (the dashed curve accounts for 
the electrostatic field a/i); b) the proton energy; c) the energy dispersion; d) the 
proton beam angular aperture. For nj = 6, wl = 10 (blue, square), n; = 12, 
wi^ = 10 (green, triangle) and rij = 6, ui^ = 20 A^ (black, circle). In panel 
b) and c), dashed curves account for predictions from the linPA model. 

equations of motion of the light ions in the aecelerating fields and atot , respectively. 
From this we obtain the traveled distance xi and the light ion layer thickness 8xi at 
the end of the linPA stage: 

XI - \/t^r+tLpA-tr, (3) 

Sxi ^di + \jtl/{l + ai/ahf +tl^pj^ - y^t^ + fLp^ - {ai/ah) tr (4) 

In what follows, we shall show that light ions exhibit two distinct long-time behaviors 
depending on whether or not they are relativistic at the end of the linPA stage. 

(i) As discussed in Sec. [21 for effective laser powers below P^'^ (which also 
corresponds to tunPA ^ tr), light ions are non- relativistic at the end of the linPA 
stage. From Eqs. ^ and ([4]), we obtain that such ions have propagated on a distance 

~ ^linPA/i'^'^r) during this first stage, while the light ion layer thickness has 
increased to Sx'l'^^ ^ di + {ai/ah)tfj^^pj^/{2tr). For most parameters of interest, this 

(c) 

thickness dxj^ remains much smaller than the transverse width wi, so that the light 

ion layer conserves its pancake-like shape. Moreover, because xj'^'' <C tunPA ^ 
light ions remains close to the heavy ion layer at the end of the non-relativistic 
linPA stage. As a consequence, they will gain further energy in the electrostatic 
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field due to heavy ions while the effect of the electron cloud on ion acceleration for 
t > tiinPA can be neglected. The acceleration process then becomes similar to DCE, 
with the difference that particles have a finite initial kinetic energy. Considering that 
both light and heavy ion layers retain their disk-like shape and that the distance 
between both layers (5a;|'^'* ^ tunPA remains small compared to the transverse 
width wj^, the accelerating field at the rear side of the heavy ion layer behaves as 
o-dce{x) = (a/i/2) (1 — 2x/wj_). In Ref. [H], the authors assume that the field 
inhomogeneity Aah ~ —ahd^i^^/w±^ through the second layer is responsible for both, 
light ions bunching and energy dispersion. However, they do not account for the 
self-consistent field ai in the second layer that, under conditions of interest in our 
work, exceeds the field inhomogeneity Aa/^ . We thus find that the field inhomogeneity 
weakly affects the long-time {t 3> tunPA) energy dispersion of light ions, which mainly 
originates from electrostatic repulsion between the protons, as previously observed in 
2D simulations (Sec. 14. ip . 

The energy gain of light ions during the DCE stage can be estimated from their 
potential energy in the heavy ion electrostatic field, while their energy dispersion is 
obtained from their potential energy in the self-consistent field a;: 



The first term in the right-hand-side of Eq. (O is similar to estimates obtained in 
Ref. [H]. The second term accounts for the fact that, in contrast to what occurs in 
"classical" DCE, part of the potential energy of light ions has already been transformed 
into kinetic energy in the linPA stage. Obviously, because light ions have not traveled 
far from the heavy ions layer during the non-relativistic linPA stage under current 
conditions, the DCE stage provides the main contribution to the final energy of 
accelerated ions. Finally, the light ion energy at < 3> tunPA is simply ^ ZiahW±/4: 
and it scales as the square-root of the laser power. Moreover, the relative energy 
dispersion ^ wiai/{w±ah) is obtained from Eqs. ([5]) and ([6|). It can be kept to a 
rather low level assuming that both the electrostatic field ai in the light ion layer and 
the transverse width wi of this layer are small compared to the accelerating field 
and the transverse width wa_ , respectively. 

In addition, we want to point out that previous considerations on 3D electron 
behavior on long times (Sec. l3.2( ) suggest that the ion beam leaves the target as a non- 
neutral bunch. The effect of electrons on long time scales, which can be disastrous in 
terms of energy dispersion, can thus be neglected: Coulomb repulsion of the light ions 
is the main source of energy dispersion. 

(ii) For laser powers above Pj^\ light ions gain relativistic energies in the 
linPA stage and the acceleration process beyond tunPA is modified. Indeed, taking 
tiinPA ^ tr in Eqs. ^ and we obtains that light ions have traveled over 
a distance a;|"'' ~ tunPA ~ tr in the linPA, while their thickness as increased 
to fe^"'^'' ^ di + {ai/ah)tr. Once more, for characteristic parameters of interest, 
^ wi, and the light ion bunch keeps its pancake-shape. However, considering 
that electrons have propagated on a distance Xe ~ tunPA, light ions should stay close 
to the electron cloud at the end of the finPA stage {x^ — ccl"''^ ^ t^ ^ w±). Therefore, 
in a second stage, light ions may gain further energy in the field due to the electron 
cloud, which is continuously pushed by the laser pulse. 

This scenario might prove very interesting for generating high-quality relativistic 



£dce ^ Zi ahWj_/4: ~ EunPA/'i. , 
AS DCE Ziaiwi/4. 



(5) 
(6) 



Ion beams by irradiating a nano- structured target with a PW laser pulse 



16 



ion beams. Nevertheless, a full 3D modeling is required to quantitatively account for 
the non-trivial long-time electron behavior. This will be the subject of future studies. 

5. Discussion and conclusion 

In this paper, we propose a novel ion acceleration mechanism allowing for an excellent 
control of all ion beam properties. This mechanism is based on complete removal 
of electrons from an ultra-thin double layer target by radiation pressure. The 
consequences of this mechanism in terms of ion acceleration are discussed here for 
the first time. Numerical simulations show that complete expulsion of electrons from 
the target can be achieved above a laser intensity threshold that depends on the areal 
charge of the irradiated target as well as on the pulse profile and polarization. For 
a nano-scale target irradiated by a CP laser pulse, laser field amplitudes ol > 50, 
typically, are required. 

We discuss acceleration of light ions of a double-layer target, in a regime where 
a capacitor-like, quasi-homogeneous, electrostatic field is built up between the ion 
and electron layers. The resulting acceleration process is referred to as linPA. It may 
accelerate light ions up to high energies over a few optical cycles only, with an energy 
dispersion that can be controlled by the target properties (areal charge and geometry) . 
Generation of relativistic light ions during this stage is also expected when using an 

(r) 

effective laser power larger than , proportional to the squared inverse charge-over- 
mass ratio. For protons, this threshold power is of a few petawatt. 

Beyond this transient linPA stage, the light ion acceleration is determined by 
multi-dimensional effects. Analytical considerations suggest that two scenarii should 
be considered depending on whether or not light ions have acquired relativistic 
velocities at the end of the linPA stage. While low energy ion generation mainly 
follows from DCE, relativistic ions are mainly accelerated in the linPA stage and may 
gain further energy in the electron cloud field afterwards. In between, modeling from 
the linPA should provide good estimates for the ion source properties. The energy 
dispersion of the resulting ion source follows mainly from Coulomb self-repulsion of 
the light ions. Energy dispersion is shown to depend mainly on the ratio of the areal 
charge in the two layers and on the target geometry. An excellent control of the energy, 
energy dispersion as well as total charge of the light ion beam can thus be achieved by 
choosing appropriate target properties. This approach thus appears to be extremely 
attractive for applications such as proton therapy, where high-quality ion beams are 
required. 

Multi-dimensional effects are partly accounted for in our theoretical model, 
but, at the present stage, we can only give rough estimates on expected ion beam 
properties. Moreover, we do not discuss the complex behavior of electrons when 
they are strongly accelerated by the laser pulse. In a realistic 3D configuration, 
inhomogeneities in the laser intensity distribution will definitely affect this process. 
As already observed by other authors [Ml [151 EQIj Rayleigh- Taylor- like or Weibel 
instabilities may occur. Another mechanism likely to influence the electron behavior 
and in turn ion acceleration is the radiation friction experienced by ultra-relativistic 
electrons in the presence of a high-amplitude laser field. All these effects on ion 
acceleration are beyond the scope of this work and are left for future investigations. 

The drawbacks of the proposed light ion acceleration scheme are similar to those 
of other methods based on the irradiation of a thin double-layer target by high power 
lasers. First, maintaining the integrity of such thin targets requires ultra-high contrast 
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laser pulses (in excess of 10^^). Also, the conversion efficiency from the laser to light ion 
is rather small and most of the absorbed laser energy goes to electrons and heavy ions. 
Moreover, techniques to separate the high-quality light ion beam for other particles 
such as heavy ions and electrons must be considered. While this is greatly simplified 
by the fact that all particles have different charge-ovcr-mass ratio, standard deflection 
and shielding techniques may drastically enhance the size of the ion source. 

Experimental exploration of this acceleration mechanism requires pctawatt lasers 
with well-controlled temporal profile and the design of nano-structured targets, such 
as diamond-like carbon foils. Fast developments in the technology of high-power 
lasers j33) and ultra-thin target fabrication [34] may make generation of high-quality, 
well-controlled, light ion beams feasible in a near future. Finally, let us note that the 
possibility to obtain high-energy ions on very short time scales (a few optical cycles) 
makes our mechanism particularly interesting for investigation at high-power, modest 
energy, laser facilities such as ELI [55] . 
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